Momentum Distribution Curves in the Superconducting State 
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We demonstrate that the E vs. k dispersion in the superconducting state extracted from momentum 
distribution curves differs qualitatively from the traditional dispersion extracted from energy distribution 
curves. This occurs because of a combination of many-body effects and the presence of an energy gap, along 
with the associated coherence factors. Analysis of such MDC dispersions can give important information on 
the microscopies of high temperature superconductors. 
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Traditionally, practitioners in angle resolved photoe- 
mission spectroscopy (ARPES) have analyzed data at 
fixed momentum as a function of binding energy, so- 
called energy distribution curves (EDCs). Recent ad- 
vances in analyzer technology have allowed the probing 



This has led to the realization that additional informa- 
tion can be obtained by analyzing data at fixed binding 
energy as a function of momentum, so-called momentum 
distribution curves. Such MDCs have been used in the 
high temperature cuprate superconductors for a variety 
of purposes, including the testing of the marginal Fermi 
liquid hypothesisa'cl, and the elucidation of a dispersion 
kink along the nodal directionci, the origin of which is 
currently being debatedETQ. 

Analysis of MDCs in the normal state, or in the super- 
conducting state along the nodal direction, is relative 
straightforward because of the absence of an energy ga; 
As we demonstrate in this paper, qualitative changes oc- 
cur in the MDCs due to the energy gap. By analyzing 
MDC dispersions, one can gain important information on 
many-body effects in the superconducting state. 

The data reported in this paper were obtained at the 
Synchrotron Radiation Center, Wisconsin, using a Sci- 
enta SESL200 analyzer, and were previously used in ear- 
lier workD. A photon energy of 22 eV was employed with 
the optimal doped (T C =90K) Bi 2 Sr 2 CaCu20 8+5 (Bi2212) 
sample in a Y — M polarization geometry. The chemical 
potential was determined from a polycrystalline Au sam- 
ple in electrical contact with the Bi2212 sample. 

As previously reportedfl, we have taken data for a num- 
ber of momentum cuts in the Brillouin zone both in the 
normal and superconducting states. In this paper, we 
will concentrate our attention on a particular momen- 
tum cut intermediate between the (tt, tt) direction where 
the superconducting gap vanishes and the (tt, 0) region 
where the superconducting gap is maximal. The reason 
for avoiding (tt, 0) is that a combination of matrix ele- 
ment effects, superstructure images, and the pseudogap, 
complicate the interpretation of MDCs in this region of 
the zoneO (for the chosen cut, these complications are not 



present). MDC and EDC dispersions were obtained from 
the maxima of the respective curves. 
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FIG. 1. (a) MDC dispersion in the superconducting state 
(SC, T=40K) versus that in the normal state (NS, T=140K). 
(b) MDC versus EDC dispersion in the superconducting state. 
k y is in units of it /a. For this momentum cut, k,;=0.597r/a. 

In Fig. la, MDC dispersions are shown in both the 
normal and superconducting states. The normal state 
dispersion is roughly linear in k in the energy range of 
interest. In the range of 20-60 meV, the superconducting 
dispersion is also linear, but with a slope approximately 
half that of the normal state, as noted earliera. This 
implies an additional many-body renormalization of the 
superconducting state dispersion relative to that in the 
normal state. Another effect of this renormalization can 
be seen at binding energies higher than 60 meV, where 
the dispersion goes almost vertical before recovering back 
to the normal state dispersion. 

To understand this effect in greater detail, we compare 
in Fig. lb the dispersions in the superconducting state 
obtained from MDCs and EDCs. As noted in an ear- 
lier paperQ, the EDC dispersions contain two branches, 
a lower binding energy "quasiparticle" branch, and a 
higher binding energy branch (known as the "hump" in 
the (tt, 0) region). We see, then, that the vertical part of 
the MDC dispersion corresponds to a crossover between 
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the low energy and high energy EDC branches. These 
effects. Me typical of electrons interacting with a bosonic 
modeErli9, and the mode in the current case, has been 

Mm and a 



identified as a spin exciton by some authors 
phonon by othersQ. As can be seen from Fig. lb and 
also noted above (Fig. la), the renormalization is an ad- 
ditional effect associated with the superconducting state, 
and thus unlikely to be due to a phonon. Moreover, in 
Fig. 2, we show the MDC and EDC dispersions calculated 
from the spin exciton modeO. We see that this calcula- 
tion gives a good description of the experimental data of 
Fig. lb, and demonstrates the pronounced effect of mo- 
mentum dependent many-body interactions on the shape 
of the MDC and EDC dispersions. This calculation also 
reproduces the difference between the MDC and EDC 
dispersions in the "linear" regime (20-60 meV), which as 
we demonstrate below, does not occur in a BCS model. 
The difference in the spin exciton case is associated with 
the mode energy scale. 
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FIG. 2. MDC and EDC dispersions in the superconduct- 
ing state from the spin exciton model of Ref. |l| (r=10 meV, 
mode energy fi= 39 meV, coupling constant g=0.65 eV, max- 
imum gap Ao=46 meV). The energy resolution is cr—1 meV 
in (a) and 7 meV in (b). In (a), the weak kink in the EDC 
quasiparticle branch marks the mode energy, which is washed 
out in (b) due to resolution. The higher energy structure as- 
sociated with the S shape in the MDC dispersion is due to the 
strong frequency dependence of the self-energy around Ao+fi. 

For the remainder of the paper, we concentrate on the 
low binding energy range, where the data are character- 
ized by a renormalized "quasiparticle" branch, and a sim- 
pler analysis of the data is possible. Returning to Fig. la, 
we note that for binding energies lower than 20 meV, the 
MDC dispersion in the superconducting state goes al- 
most vertical, and at zero energy is close to the normal 
state Fermi momentum. By looking at Fig. lb, where 
the MDC and EDC dispersions are compared, we notice 
that the upturn in the MDC dispersion corresponds to 
entering the subgap region identified from the EDC dis- 
persion. 

Unfortunately, a model independent analysis of the 
data is somewhat impractical, as noted in passing in an 



earlier paperi. This can be easily seen by a quick look 
at BCS theory. In this theory, the effect of supercon- 
ductivity can be represented by a self energy of the form 
A£/ (uj + £k + i0 + ) where Ak is the superconducting en- 
ergy gap and ek the normal state dispersion. Note that 
this self-energy has a non-trivial dependence on momen- 
tum (that is, it is not linear in k), and thus invalidates a 
simple Lorentzian analysis of the MDCs. 

We have thus looked at a simple model to describe 
the MDC superconducting dispersion in the low energy 
range. Simple BCS theory will not work, since by def- 
inition it has no solution in the subgap energy range. 
The simplest generalization is to use a self-energy of 
the formEl -iTx + A£/(w + e k + iT ). The case whexe. 
Ti = Tq = r is just a broadened version of BCS theoryO 
and we find that it gives a good account of the data. 

For Ak, we assume a d-wave energy gap of the form 
Ao(cos(k x a) — cos(fc !/ a))/2 where Ao is fit by the energy 
of the quasiparticle peak in the EDC at the Fermi mo- 
mentum. We find that £k is consistent with our earlier 
tight binding fit to normal state datalij if a scaling factor, 
z, is introduced to account for the additional many-body 
renormalization of the superconducting state dispersion 
relative to the normal state one discussed above (for the 
momentum cut considered here, z = 0.61). The peak 
energy, 24 meV, sets Ao to be 46 meV. At this stage, 
we will assume all the broadening is due to T, which is 
obtained by fitting the top of the EDC peak at the Fermi 
momentum (giving 15 meV). The effect of energy and 
momentum resolution broadening will be treated later. 
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FIG. 3. Theoretical MDC dispersion in the superconduct- 
ing state, (a) Curves correspond to Eq. 1 (BCS), second term 
in Eq. 1 only (BCSv), and Eq. 2 (phen). The circles are the 
BCS energy dispersion, uj = — Ek- The parameters are listed 
in the text (r=15 meV). (b) Results from Eq. 1 as a function 
of the broadening parameter, I\ 

In Fig. 3a, we show our theoretical MDC dispersion, 
and compare it to some alternate theoretical curves to be 
discussed below. To appreciate these results, we remind 
the reader that the broadened BCS spectral function can 
be written asu 
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where Uk and are the BCS coherence factors, and E^ = 
\/e k _ + the BCS quasiparticle energies. The simplest 
MDC is for u> = 0. In this case, the right hand side simply 
reduces to T/(r 2 + Ignoring the weak variation of 
Ak with k, one has a peak centered at £k = 0, i.e., at 
the Fermi momentum, kp. An important point is that 
for this case, the coherence factors drop out, and it is for 
this reason the peak is at kp . 

Now consider ui < (occupied states), but within the 
subgap region. If it were not for the coherence factors, 
the MDC peak would still be centered at kp, i.e., the 
dispersion would be vertical. But, the coherence factors 
skew the peak to be centered at k < kp. This trend be- 
comes more pronounced as T increases, as can be seen in 
Fig. 3b, where the MDC dispersion increasingly resem- 
bles the normal state. 

As an exercise, we show the MDC dispersion in Fig. 3a, 
but where only the second term in Eq. 1 is included. In 
this case, the MDC dispersion becomes vertical in the 
subgap region, but with an uj — value significantly dis- 
placed from kp , the displacement being due to the skew- 
ing caused by the k dependence of v-^. Physically, this 
behavior could occur if the T value for the first term in 
Eq. 1 was signficantly smaller than that for the second 
term. We note that several microscopic .theories for the 
cupratcs do in fact predict this behavioiO. That is, the 
MDC dispersion is not only a sensitive test of the coher- 
ence factors, but also particle-hole symmetry of the self- 
energy as well, even when looking at just the occupied 
states. Moreover, as we discuss below, energy resolution 
will cause this skewing to occur as well. 

Another curve is shown in Fig. 3a, and that is where Tq 
is set to 0. ThiSpSelf- energy is essentially the one used in 
our earlier workEJ in the superconducting state, and cor- 
responds to having no broadening in the BCS part of the 
self-energy. Although this model gives a good descrip- 
tion of the low energy part of the EDC at kp, it gives 
an erroneous MDC dispersion. This can be understood 
from the spectral function of this model 



irA(k, uj) 



r 2 + (w + £ k ) 2 [(c^ - E k )/(tu + e k )] 2 



(2) 



It has been factored in such a way to emphasize the u < 
branch. Note this is not of the form of Eq. 1. In particu- 
lar, the "coherence" factors are oj dependent. At u> = 0, 
the MDC is zero at kp, that is the MDC has a minimum, 
rather than a maximum as in Eq. 1. Because of this, the 
MDC dispersion has a very strange behavior in the sub- 
gap region. Moreover, for uj < —A, the MDC dispersion 
becomes similar to the EDC dispersion, reflecting the dis- 
persion of £"k- This is because of the peculiar "coherence" 
factors, which lead to a sharp second peak in the MDC 
for k > kp corresponding to the particle-hole image of 
the dispersion for positive u>, as shown in Fig. 4a. This 
subsidiary peak is strongly reduced in Eq. 1 (consistent 



with experiment), as can also be seen in Fig. 4a. To test 
this further, we have done two dimensional uj — k inten- 
sity plots, and find that Eq. 1 gives an intensity profile 
similar to experiment. This is in contrast with results ob- 
tained from Eq. 2, reminiscent of pure BCS theory (i.e., 
zero broadening), where a pronounced intensity is seen 
for k > kp, reflecting the "backbending" expected from 
the -Ek dispersion. This is not seen in experiment, and 
reiterates our point that the MDCs and intensity profiles 
are quite sensitive to self-energy effects and coherence 
factors. 
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FIG. 4. (a) Theoretical MDC in the superconducting state 
from Eq. 1 (BCS) and Eq. 2 (phen), compared to experiment 
(open circles), ui— -30 meV, and T=15 meV. (b) Same results, 
but including energy resolution (<r=7 meV, T=10 meV). 

We now discuss the effects of momentum and energy 
resolution. We have verified by calculation that the small 
momentum window of the Scienta analyzer (a rectangle 
of dimensions 0.0l7r/a along the cut, and 0.027r/a trans- 
verse to the cut) has no effect on the results. This is not 
true for the energy resolution. The latter can be deter- 
mined by fitting the leading edge of the Au spectrum, 
which for the present data gives a Gaussian a of 7 meV 
(FWHM of 16-17 meV). As an initial exercise, let us as- 
sume that all the broadening is due to energy resolution 
(for 30 meV FWHM, this would yield a a of 12.8 meV). 
Then, in the BCS case with the spectral function as delta 
functions, the MDCs are very easy to determine. For 
an energy gap much larger than temperature (satisfied 
here, since the EDC peak energy at the Fermi momen- 
tum is 24 meV, and the temperature is 40K), the Uk term 
drops out (f(Eic) is essentially zero). Doing the energy 
resolution convolution, the ARPES intensity is simply 
v 2^ e -(u+E k ) 12a _ since the intensity is totally controled 
by the ^k term, the MDC dispersion will obviously be 
skewed, as demonstrated in Fig. 5a, where this result is 
compared to the previous broadened BCS case of Fig. 3a. 
In particular, the MDC dispersion at zero energy yields 
a momentum value which is slightly displaced from the 
normal state Fermi momentum, and the MDC disper- 
sion in the subgap region is more vertical than in the T 
broadened case. 
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FIG. 5. (a) The effect of energy resolution on the MDC 
dispersion from Eq. 1. The curves correspond to T=15 meV, 
cr=0 (r), r=0, a=12.8 meV (<x), and r=10 meV, a=7 meV 
(both), (b) Comparison of theoretical (th) and experimental 
(exp) MDC and EDC dispersions. For theory, r=10 meV, 
cr=7 meV. 

Also shown in Fig. 5a is the more realistic case where 
both effects arc incorporated. Given the actual experi- 
mental a of 7 meV, a T of 10 meV is necessary to re- 
produce the 30 meV FWHM of the EDC peak. As ex- 
pected, the combined result is intermediate between the 
two limiting cases. Moreover, the inclusion of energy res- 
olution lessens the difference between the MDC profiles 
of the two self-energy models presented in Fig. 4a, as il- 
lustrated in Fig. 4b. An advantage of the model of Eq. 2 
is that it can account for the extra experimental weight 
on the trailing (unoccupied) edge of the MDC peak not 
present in the broadened BCS model. In fact, by looking 
at experimental MDCs at higher binding energies, a weak 
shoulder does develop on the trailing edge, correspond- 
ing to the expected particle-hole image discussed earlier. 
This image should become better defined with improved 
resolution and statistics. 

In fact, it is somewhat remarkable that the simple form 
of Eq. 1 does such a good job in describing the data. In 
Fig. 5b, we compare the MDC and EDC dispersions of 
Fig. lb to our calculation of Fig. 5a. Confining ourselves 
to the "quasiparticle" branch, the agreement of exper- 
iment and theory is quite good. But, the discrepancy 
between the MDC and EDC dispersions in the "linear" 
regime (20-60 meV) is not reproduced by BCS theory, 
though it can be accounted for in the spin exciton model 
as demonstrated in Fig. 2. 

In conclusion, we find that the MDCs and resulting dis- 
persions are non-trivial in the superconducting state, and 
give important information on the electron self-energy 
and coherence factors. We feel that a more detailed study 
of MDCs, both in the superconducting and pseudogap 
phases, will give important insights into the microscop- 
ies of high temperature cuprate superconductors. 
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